Transitions involving conical magnetic phases in a model with bilinear and 

biquadratic interactions. 
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In a previous work a model was proposed for the phase transitions of crystals with localized mag- 
j^Tjj netic moments which at low temperature have a "conical" arrangement that at higher T transforms 

^ , into a more symmetrical structure (depending on the compound) before becoming totally disordered. 

The model assumes bilinear and biquadratic interactions between magnetic moments up to the fifth 
neighbours, and for any given T the structure with the least free energy is obtained by a mean-field 
approximation (MFA). The interaction constants are derived from ab initio energy calculations. 

In this work we improve upon that model modifying the MFA in such a way that a continuous 
(instead of discontinuous) spectrum of excited states is available to the system. In the previous work, 
which dealt with LaMn2Ge2 and LaMn2Si2, we found that transitions to different structures can be 
. obtained for increasing T, in good qualitative agreement with experiment. The critical temperatures, 

' however, were exaggerately high. With the new MFA we obtain essentially the same behaviour 

concerning the phase transitions, and critical temperatures much closer to the experimental ones. 

e 
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Localized magnetic moments in several compounds of the family RT2X2 (R = rare earth, T = transition metal, 
q \ X = Ge, Si) are found to have "conical" ground state (GS) configurations (see eq.l) and to undergo transitions to 
configurations with higher symmetry (helical, canted, collinear) as temperature T increases [l]. Among the many 
works on this subject we can cite refs.2-6 (see also references therein and in ref.l). In this work we consider LaMn2Ge2 
and LaMn 2 Si2 In such compounds only the Mn atoms bear a localized magnetic moment, so that we consider 

only the Mn sublattice. We take it as a base-centered tetragonal lattice, whose lattice constants are given, in terms 
of the crystallographic constants a and c, by at = a and cf = c/2, and the lattice sites by Rhki = haie x + kafe y + lcfe z , 
where I is an integer and j, k are both integers or half-integers. In the following we will use subindices i and j as short 
notations for sets (h, k, I). We work with the total atomic angular momentum (in units of %), which we call spin and 
denote S, as in the spin Hamiltonian formalism. The corresponding magnetic moment is jl = (J,bS ( is Bohr's 
, magneton). 

The magnetic structure can be characterized by a pair of polar angles 9 and a, such that at any site Rhki = 
• • . ha'e x + ka'e y + lc'e z the average value of the local spin is given by 
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I. INTRODUCTION 



Shki = S[sin(9)cos(la + £hk)e*x + sin(9)sin(la + £,hk)e y + cos(9)e z ], (1) 

where = or n for h, k integers or half-integers respectively and the bars indicate the statistical averages (see also 
eqs. 22-24). We can also say that the polar angles (#,, <f>i) characterizing the direction of the average spin at a site 
i = (h, k, I) are given by 9i = 9 and </>j = (la + £hk)- For general values of 9 and a (i.e. values different from 0, ^7r, it 
and |7r) this gives what is customarily called a conical structure (even though perhaps the term "conical" tout court 
should be limited to the case £hk = V h, k). Particular cases of the conical structure are the helical (9 — \tt), the 
canted ( a = 0, \tt, it or §7r") and the collinear (all the spins are parallel or anti-parallel to each other) structures. 
In the following we will call "general conical structure" a conical structure in which 9 and a can take any values, 
including these particular values. 

As T increases the spins get progressively disordered (so that their average values decrease), while angles 9 and a 
vary, until eventually taking values corresponding to some of the above-mentioned "particular cases" so that a more 
symmetrical structure is obtained. In a previous work[7| a model was proposed assuming that the localized magnetic 
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moments are subjected to bilinear and biquadratic interactions up to fifth neighbours (and possibly to an anisotropy 
field). The interaction constants were obtained from ab initio calculations of the energies of different structures. 

For any given T the structure was obtained by finding, in a mean- field approximation (MFA), the 9 and a angles, the 
modulus of the average spins, and the average of to 2 (to = total angular momentum quantum number) for which the 
free energy was minimal. This model could qualitatively account in a reasonable way for the mentioned behaviours, 
but it gave exaggerately high values for the transition temperatures T c 's between the different configurations (see 
sections IV, V). In this paper we modify the MFA in such a way that reasonable values of the T c 's are obtained. 

II. THE MODEL 

A common approach to a phenomenological study of a crystal with localized magnetic moments is the employ of 
the bilinear Heisenberg Hamiltonian: 

Hj = ^ JijSi ■ Sj (2) 

(the prime in the summation means i ^ j). 

The ground state for such a Hamiltonian is a general type of helical structure Q , so that in order to obtain a conical 
structure some other interaction must be added to it. The addition of a biquadratic term to the Heisenberg bilinear 
interaction was shown to explain another type of phase transitions in a RT2X2 compound, namely U Ni2Si2 @. In 
ref.7 it was found that a Hamiltonian with bilinear and biquadratic interactions both up to the fifth neighbouring 
sites gives an acceptable fitting of ab initio total energy values for several configurations in which the moduli of the 
spins are all equal to S. The anisotropy field was not required, but we included the quadratic term in order to study 
its possible effects (a cos 2 (#) term does anyway appear in the bilinear interaction). 

Accordingly, we assume the Hamiltonian 

H = Hj + H B + H D (3) 

with 

h b ^-\y,' b ^-^ 2 ( 4 ) 
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and 



H D = D 2 J2sf. z (5) 

i 

where the Jij and are non-zero for the five nearest sets of neighbours (which, for the reference site (0,0,0), are: 
(±±±±,0); (±1,0,0) and (0, ±1,0); (0,0, ±1); (±1, ±1,0); (±±,±±,±1) 

Let us denote | rrii > the eigenstate of a spin at site i with eigenvalue rrii in the direction given by the polar angles 
(Ouch). 

The interaction energy between the spins at sites i and j (i.e. 
< m,i |< m,j | (Hj + Hb) | >| rrij >) is 

UjB;ij(mi,mj) = Uj- l j{m t ,m J ) + V (m l; m^), (6) 

with 

Uj ;ij (mi,mj) = -yijiii.tn ,i>,, (7) 

and 



-P(S*+ 1) - 3to2]P(S +1) - 3m 2 j](pl + 1) + S{S + \){mj + m)) - 3m 2 m 2 - -m.m^ 



(8) 



3 



where pij is the cosine of the angle between the spins at sites i and j. The matrix element of Hb was calculated 
by rotating the spin operators [Tol ] . 

In a general conical structure all sites are energetically equivalent to each other, so that neglecting boundary effects 
the total energy is equal to NU S , where N is the total number of sites and U s is the energy per site. As reference site 
we take (h = k = I = 0), also denoted i = 0. We have 

N-l 

Us = ]T UjaMmw) + -D 2 [[5(5 + 1) - ml] + [3m 2 - 5(5 + 1)] cos 2 (0)] . (9) 
i=i 

We shall denote respectively J„ and B v the constants Jy and Bij for sites i,j that are z/ h neighbours of each other. 
In ref.7 eq.(9) with mo = rrij — 5 was used to determine the parameters {J V ,B V ,D 2 of the Hamiltonian. The 
energy dependence on 9 and a was written as 

U T=Q (9,a) =Xx + A 2 cos 2 9 + A 3 sin 2 6 cos a + X 4 cos 4 + A 5 (cos 2 + sin 2 9 cos a) 2 + X 6 (cos 2 9 - sin 2 9 cos a) 2 (10) 

where the A„'s are constants. 

In principle the A„'s can be obtained by fitting the results of ab initio total energy calculations for several config- 
urations with m = rrij = 5. Such energies were obtained by using the FLEUR code, as described in ref.7. If this 
fitting is done for six configurations, i.e. for six pairs (9, a), a linear system is obtained, which can be solved for the 
X n . We did this for several sets of six configurations, so as to check the consistency of the formalism (i.e, that in 
each case the same set of X n is obtained with a reasonable accuracy). We found that this consistency is achieved if 
these configurations are neither too close to each other nor too far from the GS. In the former case the determinants 
appearing in the numerical solution of the system are nearly zero (they are sums of terms which nearly cancel out 
each other) ; in the latter case it is possible that if one forces the electrons associated with the magnetism to occupy 
high-energy states, one cannot neglect the additional energy due to changes of the other electronic states. In both 
cases the error is large. 

The A n 's to be adopted in this work for LaMii2Ge2 and LaMn2Si2 are the averages of the A„'s of the selfconsistent 
sets. 

Having obtained the X n 's, we must face the problem of determining the eleven parameters appearing in the 
Hamiltonian. 

The part of the magnetic energy not depending on 9 and a (i.e. X{) is not the same as the analogous part of the 
total energy (as a matter of fact, it is several orders of magnitude smaller), so that it cannot be singled out in the ab 
initio values. Therefore, X\ is of no use for the determination of the parameters. As a consequence, J 2 , J4, B 2 and B^ 
cannot be determined because they multiply terms with ^02 or po4 5 which in a conical structure are constant, being 
equal to 1. In order to express the other seven parameters in terms of the five constants X n (n=2,. . .,6) we compare 
eq.(10) with the energy per site (eq.9 with (6-8)) for m, = m 3 = 5 V In doing so we must take into account that 
any site has four first, second and third neighbours, two fourth neighbours and eight fifth neighbours. We obtain the 
system 

X 2 = -i5 2 [4(2J x - BO + (2J 3 - S 3 ) + 4(2J 5 - B 5 )} + 85 2 (5 - ^) 2 B, + 5(5 - i)D a 
A 3 = -^S 2 (2J 3 -B 3 -8J 5 +AB 5 ) 

A 4 = -85 2 (5- i) 2 Sx] (11) 
A 5 = -5 2 (5-i) 2 B 3 
A 6 = -45 2 (5-i) 2 5 5 

From here we obtain directly B±, B3 and B§ , after which we are left with two equations for the unknowns J\, J3, 
J5 and D 2 , two of which must be chosen arbitrarily (as well as J 2 , J4, B 2 and B4). Concerning D 2 , we shall give it 
a small value, while for the other three unknowns we think that the less arbitrary procedure is that of choosing J5 
in such a way that it and the values of J\ and J3 calculated in terms of it and of D 2 have a smooth variation with 
the interatomic distance, and are, in principle, about an order of magnitude greater than the corresponding Bj. An 
analogous criterion shall be followed for J 2 , J4, B 2 and -B4. 



Setting ±S 2 (S* - ±) 2 = C 2 we have 



Bi = 


16 ' 


(12) 


Bi — 


2 5 ' 2 


(13) 


B 5 = 


-\x 6 /c 2 


(14) 


J 3 = 


-\x 3 /S 2 + ^B, + J 5 -2B 5 


(15) 


Ji = 


\-X 2 /S 2 + [8(5 - i) 2 + 2]Bi - J 3 + ^3 - 4J 5 + 2S 5 + [(S - \)/S]D 2 


(16) 



III. CALCULATION OF THE FREE ENERGY 



In this work we must find the magnetic structure of the LaMn2X2 compounds (which we know from experiment 
to be a general conical structure) at any temperature T. Therefore we must determine which is for the different sets 
of values of the interaction constants the conical structure that minimizes the free energy. In a general state of the 
system the spin at each site i is in a state | m, >, so the system is characterized by the 3iV quantities TOj, 9i, 4>i 
(i = l,...,N). 

For T > the state of the system will be that superposition of the above-defined states which minimizes the free 
energy F = — ksT log Z, i.e. which maximizes the partition function. The latter is given by 



Z= / dn o- / dfijv-iexp[-/?l/(m i ,0 i ,&)] ) (17) 

mo,.-) m JV — 1=— S 

where (3 == 1/kT, dili = sin 8 idOidfa and each integration is carried out over all the space directions. 

Taking advantage of the fact that all the sites arc physically equivalent to each other we define the partition function 
per site Z s through Z = Z^ . 

We calculate Z s for the reference site i = 0. To do this we must make a MFA, i.e we must approximate the energies 
of the states at that site in the field of the spins at the other sites (which will be referred to as "source sites") by 
replacing each of the latter by some average state. This is done as follows. The rrij , 6j and otj at each source site j 
are taken according to eq.(l), where S, 9 and a are unknowns. Furthermore, as the contribution of a j site to the 
energy of the reference spin depends also (due to the biquadratic interaction) on m 2 we must assign to each source 

site a fourth variable, to 2 . We have 



s , 

Z S MFA(0,a, S, to 2 ) = Y m dn exp[~j3U(m ,6 ,(l)o;0,a,S,m 2 )]. (18) 



mo — — S 



The energy in the integrand contains the contributions of all the source sites, and consists of an interaction term 
and an anisotropy term (eqs.7,8,9). The latter, being a local term, requires no approximation. 
For the interaction energy associated with the reference site, whose exact value is 



N-l 

^UjB-fljimo^jipoj) (19) 

(where UjB;Oj( m a,'>Tij'TPnj) is given by eq.(6) with (7,8) taking into account that the spin at the i — site has the 
general orientation (9 , <fio) ), there are several possible MFA expressions, and it is not clear which one is the best in 
any particular case (see section VI). 

Now we must find, for any given T, the values of 0, a, S, and to 2 characterizing the state of minimal free energy, 
i.e. of maximal Z . For a conical structure the cxtremum conditions are 
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dZ/89 = (20) 
dZ/da = 0. (21) 

On the other hand, both S and m 2 must satisfy a self-consistence condition, i.e. they must be equal to the respective 
average values they give rise to at the reference site. S is the modulus of 

S = Z s 1 mo / <i£!o[sin6>o cos (fro e x + sin 6>o s'mcfro e y + cos 9o e z ] exp[— (3U(mo, 0q, (fro; 9, a, S, m 2 )], (22) 

TTLQ — — S 

i.e. is given by 

§ = {s; 2 + s; 2 }i (23) 

(due to the symmetry of the structure, the y component vanishes), while 
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m 2 = Z~ 1 dn m 2 exp[-f3U(m Q ,6o,cj) Q ;6,a,S,m 2 )}. (24) 

mo=—S 

So, we must solve the system of equations (20), (21), (23) and (24) [la ]. 

For the more symmetrical structures the number of equations is reduced: for the helical structure we set 9 — hix 
and drop eq.(16a); similarly, for the canted structure, we set a — n, and instead of (16a) we drop (16b). Finally, in 
the collinear structure we have 9 = ^tt, a — n, and the equation system is reduced to (14) and (15). 



IV. NUMERICAL CALCULATIONS. LAMN 2 GE 2 



The interaction constants Jij,Bij and D 2 are obtained from the X n 's used in ref.7, taking into account that the 
latter were evaluated for a cell with four sites, i.e. must be divided by 4. According to the ab initio calculations of 
ref.7 and to the experimental results of ref.2 (where the helical and FM components of the Mn magnetic moment are 
found to be approximately 2.6/xb and 1.61/Xb respectively), and assuming L = because of quenching, we take S = I. 

As we said above, only B\, £? 3 and B 5 can be derived from the ab initio energies. It turns out that these B u 's are 
negative, and their absolute values decrease smoothly for increasing interatomic distance. 

For the other constants we proceed as explained in section 2: we choose first J5 and D 2 , thereby determining Ji and 
J 3 , and then, without affecting the other constants, we choose J 2 , J4, B2 and B 4 . We can get a set of J„'s about an 
order of magnitude greater than the B^s o nly if J\ and J3 have different signs. This is a quite reasonable condition, 
being a property of the RKKY interaction [lj] 

Among the sets of parameters that we used in our calculations, one of those yielding the closer agreement with the 
experimental results is: Ji = -0.002137, J 2 = -0.0005, J 3 = 0.002792, J 4 = 0.0010, J 5 = 0.000755, 

Si = -0.000482, B 2 = -0.00023, B 3 = -0.0001555, B 4 = -0.00009, B 5 = -0.000053, D 2 = 0.0005. 

With these parameters we have at low-temperatures (LT) a conical structure with 9 ~ 60° and a ~ 114°. 

The experimental values are 9 ~ 58° and a ~ 133° (see ref.2, taking into account that therein a refers to the 
rotation of the spin at site (5, ^, 1) (instead of (0,0,1), as in this work) with respect to that at (0,0,0), and that their 
lattice constant in the z direction is c while ours is sc. 

As T increases, both 9 and a mcrease, until the helical structure is attained (with a ~ 144°). This happens at 
T c i ~ 439 K. Meanwhile, S and m 2 decrease from their LT values 1.5 and 2.25. Upon a further increase of T (keeping 
9 = 90°) a keeps increasing, reaching 180° (collinear structure) for T c2 ~ 448 K. 

Finally, both 9 and a are kept constant at 90° and 180° respectively. As T increases, S and m 2 decrease until, for 
T = T C 3 ~ 1041 K, the former gets equal to zero (paramagnetic structure). Another set of parameters obtained from 
the same X n 's is 

Ji = -0.00223, J 2 = -0.0005, J 3 = 0.00297, J 4 = 0.0010, J 5 = 0.0008, 

B l = -0.000482, B 2 = -0.00034, B 3 = -0.0001556, B 4 = -0.0001, B 5 = -0.000053, D 2 = 0.0005. 

The results are qualitatively the same as for the previous set, and the T c 's are 442.5, 455.3 and 1116 K. Notice 
that while T c \ and T c2 decrease with respect to the previous ones, T c3 increases (a behaviour we found in many other 
cases). 
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Experimentally, the transition "conical < — > helical" is indeed observed ( T c \ ~ 320 K), but from the helical 
the system goes to the PM structure ( T c2 ~ 420 K) skipping the collinear structure [<|]. However, the complete 
sequence of transitions obtained in this work is observed in a slightly different compound, namely Lai_ :E Y a ;Mn2Ge2, 
for 0.1 < x % 0.2 (see fig. 13 of ref.4). There, as in our results, the temperature range of the helical structure is 
relatively small. 

It is worthwhile to notice that in the PM phase the spins are not totally disordered (i.e. the states with m = ±i 
are less populated than those with to = ±|, so that m 2 is not equal to 1.25) until a much higher T is reached. 

V. NUMERICAL CALCULATIONS. LAMN 2 SI 2 

As for LaMn2Ge2, the interaction constants are obtained from the X n 's used in ref.7. Concerning the S of Mn, 
now both ab initio calculations 0, [HJ and experimental results give a magnetic moment near to |/iB 5 which by 
Lande formula and for a quenched orbital momentum gives a total angular momentum (denoted S in this work) equal 
to 1. In this case no conical structure is obtained for our Hamiltonian. However, ab initio studies[l2j indicate that 
this lowering of the magnetic moment is due to Si-Mn hybridization, so that, as for X = Ge, we have S = 1.5. 

Another possibility, leading to the same S, is that of having an incomplete quenching, L = 1, with the total spin 
momentum and the total angular momentum both equal to |. This would give a Lande factor g = 26/15 i.e. a 
magnetic moment /i = 2.6^, in good agreement with the above value. 

Now the set of parameters which we found to give the closest agreement with experiment is: 

Ji = -0.00145406, J 2 = -0.0003, J 3 = 0.00177942, J 4 = 0.0010, J 5 = 0.0006, 

B 1 = -0.00045278, B 2 = -0.0003, B 3 = -0.00024431, B 4 = -0.00017, B 5 = -0.00010941, D 2 = 0.000. 
The behaviour of the J v , B v in terms of v is very similar to that of the Ge compound. 

Proceeding as for the latter we have at LT a conical structure with 9 ~ 53° and a ~ 139° The agreement of these 
angles with experiment is not as good as for the Ge compound [15| . 

As T increases, both 8 and a increase, but now the former does it more slowly, so that a reaches 180° (i.e. the 
canted structure is attained) while 9 is still nearly 54°. This happens at T cl ~ 269 K. 

When T (keeps increasing (with a = 180°) 9 increases more rapidly, reaching 90° for T = T c2 ~ 393 K. This 
sequence of transitions is the one observed experimentally Q. 

Finally, in the collinear structure S and to 2 decrease with increasing T, until, for T = T c3 ~ 1061 K, the former 
gets equal to zero (paramagnetic structure). The experimental values reported in ref.5 are T c i = 45 K and T c2 = 305 
K (the T > 305 K region is not studied therein). 

VI. ABOUT THE MFA'S. 

There are many ways of making a MFA in systems like those we study here. One can approximate the exact 
expression of the energy of the spin at the reference site (i = 0) in the field of the spins at the other sites (j) in 
different ways, and thereafter replace nij and m 2 by their average values S and to 2 ) respectively. Some of such 
approximations are the following: 

- since in the expression for the total energy the i and j indices are dummy, one can replace the terms with (m 2 +m 2 ) 
in the summation over i,j by 2to 2 or 2to 2 before of choosing site i = as reference; 

- one can replace in H j each spin operator Si by Si + A j with A, = Si — Si, and neglect higher-order terms in the 
A operators Q; 

- one can do the same thing for the rmmj term in (8); 

-2 

- for the to|to 2 one can proceed analogously to the case of bilinear terms introducing the operators A2 ; ; = Si — mf. 

Having chosen the approximate expression of the energy, one can make different choices for the excited states 
available to the spin (see below). 

In this work we found that one can obtain completely different results with different MFA's. The most interesting 
example of this (and the only one for which we have an explanation!) is the spectacular decrease of the calculated 
critical temperatures when the spin at the reference site is allowed to have a set of excited states with a continuous 
energy spectrum, instead of a discontinuous one. 

In ref.7, with a certain set of parameters and a MFA using a discontinuous spectrum, we obtained for the T c 's of 
LaMn2Ge2 the values 1943 K, 2318 K and 4211 K. With the same parameters and one of the MFA's used in this 
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work (continuous spectrum) the corresponding values are 468, 473 and 977. A similar situation occurs for the Ge 
compound. This difference in the TVs can be explained as follows. 

In ref.7 the spin at site i could be in any of the eigenstates | m > with m = ±|,±| with respect to the (6, a) 
direction corresponding to the directions of the other spins. As for any given T 9 and a, as well as S and m 2 ), are 
fixed, for any set of values of these four quantities the spin at the reference site needs a finite energy in order to jump 
to an excited state. 

In the present work, as described in section III, that spin can have (in addition to the different m values) any 
orientation (#o,</>o)j i-e it can vary its orientation by an infinitesimal angle, so that it can acquire a small degree of 
desorder at a very low T. This disordering causes a decrease in S, which in its turn reduces the field at the reference 
site, making easier for the spin to jump to an excited state. In addition to this, at higher T"s the m = S — 1 state 
begins to be populated, so that also m 2 starts to decrease from its LT value S 2 . This feedback process must be what 
causes the great decrease of the calculated T c 's. 

Concerning the way of calculating the energy of the state at the reference site, in this work we tried several 
formulae, i.e. several MFA's, and ended up using three of them (that we call MFAO, MFA2 and MFA4). Let us 
consider ^\ . UjB-.ij ( m which, however, as what we need is the energy per site the summation over i is not carried 
out). . 

In MFAO we simply replace (raj and mf) by S and m 2 ) respectively. 

In MFA2 we replace J2 i ■ Bij{m 2 + m 2 )((pij 2 + 1) by j Bijm 2 ((pij 2 + 1), after which we take i as the reference 
site (i — 0). This is a usual procedure in order to have a more exact value of the total energy before making a MFA. 
In this case, in which we do not carry out the summation over i, it is not clear this to be an improvement. 

In MFA4 we use the A and A 2 operators for the terms bilinear and biquadratic in m, and m,j, after which we 
proceed as for the MFA2. 

In all cases the spin at site i is given as described in section III. 

In the calculations we carried out for the Ge compound the sequence "conical < — > helical < — > collincar < — ► PM" 
was obtained for several sets of parameters with the MFA2, but not with the MFA4 (where in most cases the system 
of equations to be solved (section IV) ceases to have a solution above a certain temperature). 

Surprisingly, for the Si compound something nearly opposite is true: the sequence "conical < — > canted < — > collinear 
< — > PM" was obtained for three sets of parameters with the MFA4 and only for one set with the MFA2 (in the latter 
case the free energy is discontinuous at T = T c2 ). 

With the MFAO we obtained the mentioned sequences for only one set of parameters for each of the compounds (in 
the case of Si, with the discontinuity at T C 2). 

VII. CONCLUSION. 

This work is concerned with the hypothesis that the great variety of magnetic structures observed in different 
RT2X2 compounds at different temperatures can be explained by assuming that the localized magnetic moments are 
subjected to bilinear and biquadratic interactions and that the interaction constants can be derived from ab initio 
calculations of the energies of different unstable structures. The following results of this work give support to this 
hypothesis: 

- starting from ab initio energies we find a set of interaction constants for which the observed structures and the 
transitions between them are predicted for LaMn 2 Si2; 

- the interaction constants we obtained by the same procedure for LaMn 2 Ge2 yield the observed structures and the 
transitions between them of a very similar compound, i.e. La .8Y .2Mn 2 Ge2; 

- for both compounds the dependence of the constants on the interatomic distance is perfectly sound for both the 
bilinear and biquadratic interactions; 

- due to the many competing interactions, the behaviour of the system is very sensitive to small changes of the 
constants, so that it can be expected that the same model can be applied to other compounds of this family, for which 
different magnetic structures have been observed. 

If our hypothesis is valid, this paper is only a first step towards a complete theory. Among the additional studies 
that should be carried out we can list the following. 

Further work is required to improve the calculation of the interaction constants. 

It is important to understand the pros and cons of each MFA, in order to adopt the most adequate one. 
It is likely that more neighbours must be included in the calculation. 

One should also take into account a possible variation of the interaction constants with temperature due to the 
change of the lattice constants, which implies a variation of the interatomic distances. 
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